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The ferromagnetic Kondo model with classical corespins is studied via unbiased Monte-Carlo sim- 
ulations. We show that with realistic parameters for the manganites and at low temperatures, the 
double-exchange mechanism does not lead to phase separation in one-dimensional chains but rather 
stabilizes individual ferromagnetic polarons. Within the ferromagnetic polaron picture, the pseu- 
dogap in the one-particle spectral function Ak(uj) can easily be explained. Ferromagnetic polarons 
also clear up a seeming failure of the double-exchange mechanism in explaining the comparable 
band widths in the ferromagnetic and paramagnetic phase. For our analysis, we extend a simplified 
model, the finite temperature uniform hopping approach (UHA), to include polarons. It can easily 
be evaluated numerically and provides a simple quantitative understanding of the physical features 
of the ferromagnetic Kondo model. 
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I. INTRODUCTION 

Manganese oxides such as Lai-^Sr^MnOa and 
Lai_2;Ca x Mn03 have been attracting considerable at- 
tention since the discovery of colossal magnetoresistance 
(CMR)ii£. These materials crystalize in the perovskitc- 
type lattice structure where the crystal field breaks the 
symmetry of the atomic wave function of the manganese 
e?-electrons. The energetically lower t 2g levels are oc- 
cupied by three localized electrons. Due to a strong 
Hund coupling their spins are aligned, forming a localized 
corespin with S = 3/2. The electron configuration of the 
Mn 3+ ions is t 2g eg, whereas for Mn 4+ ions the e g electron 
is missing. Due to a hybridization of the e g wave function 
with the oxygen 2p orbitals, the e g electrons are itinerant 
and can move from an Mn 3+ ion to a neighboring Mn 4+ 
via a bridging 2 ~. The interplay of various physical 
ingredients such as the strong Hund coupling (Jh) of 
the itinerant electrons to localized corespins, Coulomb 
correlations, and electron-phonon coupling leads to a 
rich phase diagram including antiferromagnetic insulat- 
ing, ferromagnetic metallic and charge ordered domains. 
The dynamics of the charge carriers moving in the spin 
and orbital background shows remarkable dynamical fea- 
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Since full many-body calculations for a realistic model, 
including all degrees of freedom, are not possible yet, sev- 
eral approximate studies of simplified models have been 
performed in order to unravel individual pieces of the 
rich phase diagram of the manganites. The electronic de- 
grees of freedom are generally treated by a Kondo lattice 
model which in the strong Hund coupling limit is com- 
monly referred to as the double-exchange (DE) model, a 
term introduced by Zener— . In addition, the correlation 
of the itinerant e g electrons is well described by a nearest 
neighbor (n.n.) Coulomb interaction. The on-site Hub- 
bard term merely renormalizes the already strong Hund 
coupling. For the Kondo model with quantum spins it is 
still impossible to derive rigorous numerical or analyti- 
cal results. If the S=3/2 corespins are treated classically, 



however, the model can be treated by unbiased Monte 
Carlo techniques. The impact of quantum spins on the 
electronic properties has been studied in Ref. It 
appears that quantum effects are important for (S=l/2) 
corespins or at T = 0. For finite temperature and S=3/2, 
classical spins present a reasonable approximation. 

Elaborate Monte Carlo (MC) simulations for the FM 
Kondo model with classical t 2g corespins have been per- 
formed by Dagotto et aZi 9 i 10 ' n , Yi et alH, and by Fu- 
rukawa et a /J&14. Static and dynamical properties of 
the model have been determined. These studies revealed 
features (discontinuity of the mean electron density as a 
function of the chemical potential, infinite compressibil- 
ity) which have been interpreted as signatures of phase- 
separation (PS). PS has also been reported 1 ^ from com- 
putations based on a dynamical mean field treatment of 
the DE model at T = 0. A phase diagram and critical 
exponents of the DE model have been determined with 
a Hybrid MC algorithmic. 

In the manganites, the Hund coupling Jh is much 
stronger than the kinetic energy. Consequently, config- 
urations are very unlikely in which the electronic spin 
is antiparallel to the local corespin. The present au- 
thors have proposed an effective spinless fermion (ESF) 
model— that takes effects of antiparallel spin configu- 
rations into account via virtual excitations. It has been 
demonstrated that the results of the ESF model are in ex- 
cellent agreement with those of the original Kondo model 
even for moderate values of Jh- This applies also to fea- 
tures which have been previously interpreted as signa- 
tures of PSii. Taking Coulomb interactions into account, 
PS has been argued to lead to either small 19 or large^S 
(nano-scale) clusters, which have been the basis for a 
possible though controversial 6 explanation of CMR 21 - 22 . 
Moreover, lattice distortions^ 3 * 2 ^ are believed to play a 
crucial role for the CMR effect*^ and should also be 
included in the model. 

In this paper, we present a numerical and analytical 
study of the ID ferromagnetic Kondo model with clas- 
sical corespins. We find that the correct physical inter- 
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pretation of the features which have been interpreted as 
PS in the one-dimensional model, is rather given by fer- 
romagnetic polarons, i.e. small FM-regions with one sin- 
gle trapped charge-carrier, compatible with exact diago- 
nalization results for small clusters^. This applies even 
without n.n. Coulomb repulsion invoked in Refs H9II20I . 
Energetically, there is no significant difference between 
polarons, bi-polarons or even charge accumulations in the 
PS sense. It is rather the entropy, which even near zero 
temperature clearly favors polarons. The polaron picture 
allows also a straight forward and obvious explanation of 
the pseudogap, which has been previously observed in 
the spectral density 18 i 21 i 26 i 27 . 

In a previous paper— we introduced the uniform hop- 
ping approximation (UHA), which replaces the influence 
of the random corespins on the e g electron dynamics by 
an effective uniform hopping process. Essential physical 
features of the original model could be described even 
quantitatively by UHA, while the configuration space, 
and hence the numerical effort, was reduced by several 
orders of magnitude. Besides the numerical advantage, 
UHA also allows the derivation of analytical results in 
some limi ting cases at T = 0. 

In Ref. |23 we have extended UHA to finite tempera- 
tures by allowing for thermal fluctuations of the uniform 
hopping parameter. By taking into account the density 
of corespin states, it is possible to calculate thermody- 
namic quantities of one and three-dimensional systems in 
the UHA. The reliability of this finite-temperature UHA 
has been scrutinized by a detailed comparison of the re- 
sults for various properties of the ferromagnetic Kondo 
model with unbiased MC data in ID. 

Here, we will generalize UHA to regimes, where a sin- 
gle hopping parameter is not sufficient to describe the 
physics of the FM Kondo model. Particularly near half 
filling, a typical corespin configuration shows small fer- 
romagnetic domains (polarons) immerged in an antifer- 
romagnetic background. Therefore, two different UHA- 
parameters are necessary to model the impact of the fluc- 
tuating corespins on the e g electron dynamics. 

This paper is organized as follows. In Sec.|H]the model 
Hamiltonian is presented and particularities of the MC 
simulation for the present model are outlined. The gen- 
eral discussion of ferromagnetic polarons near half filling 
is presented in Sec. lIIII Section llVl develops a generaliza- 
tion of the uniform hopping approach in order to treat 
FM polarons. Polaronic features in the spectral density 
are analyzed. The key results of the paper are summa- 
rized in Sec. [V] 



II. MODEL HAMILTONIAN AND UNBIASED 
MONTE CARLO 

In this paper, we will concentrate solely on properties 
of the itinerant e g electrons interacting with the local t2 g 
corespins. We also neglect the degeneracy of the e g or- 
bitals. The degrees of freedom of the e g electrons are 



then described by a single-orbital Kondo lattice model 28 . 
As proposed by de Gennes 2 ^, Dagotto et aZ£21 and Fu- 
rukawa-2, the ti g spins Sj are treated classically, which 
is equivalent to the limit S — > 00. The spin degrees of 
freedom (S) are thus replaced by unit vectors Si , param- 
eterized by polar and azimuthal angles 9i and 4n , respec- 
tively. The magnitude of both corespins and e g -spins is 
absorbed into the exchange couplings. 

A. Effective Spinless Fermions 

It is expedient to use the individual t2 g spin direc- 
tion Si as the local quantization axis for the spin of the 
itinerant e g electrons at the respective sites. This rep- 
resentation is particularly useful for the Jh — > 00 limit, 
but also for the projection technique, which takes into 
account virtual processes for finite Hund coupling. As 
described in Ref. ITsl the energetically unfavorable states 
with e g electrons antiparallel to the local ti g corespins 
can be integrated out. This yields the ID effective spin- 
less fermion model (ESF) 

<i,j> i,j ' <i>3> 

(1) 

The spinless fermion operators c- correspond to spin-up 
electrons (relative to the local corespin-orientation) only. 
The spin index has, therefore, been omitted. With re- 
spect to a global spin-quantization axis the ESF model 
still contains contributions from both spin-up and spin- 
down electrons. 

The first term in Eq. corresponds to the kinetic en- 
ergy in tight-binding approximation. The modified hop- 
ping integrals t"'? depend upon the t2 g corespin orien- 
tation 

= f o <f . (2) 

where the relative orientation of the ti g corespins at site 
i and j enters via 

<f(S)=cos(^/2)e^ 
<^(5)=sin(^/2)e^ ' 

These factors depend on the relative angle i?y of 
corespins Si and Sj and on some complex phases i/'ij 
and Xij- 

The second term in Eq. Q accounts for virtual hop- 
ping processes to antiparallel spin-corespin configura- 
tions and vanishes in the limit Jh — > 00. The last term 
is a small antiferromagnetic exchange of the corespins. 

It should be noted that the unitary transformation to 
the local spin quantization axis is not unique. This fact 
can be exploited to eliminate the phase factors ip+j in ID. 
Then the n.n. hopping integrals are simply given by the 
real numbers cos(-dij /2). 
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B. Grand Canonical Treatment 



We define the grand canonical partition function as 



Z = J V[S] tr c e 

[ V[S\ = TT ( / dd * sind 
J i=1 W 



2tv 



(4) 



where tr c indicates the trace over fermionic degrees of 
freedom at inverse temperature /3, N is the operator for 
the total number of e g electrons and [i stands for the 
chemical potential. Upon integrating out the fermionic 
degrees of freedom, we obtain the statistical weight of a 
corespin configuration S that can be written as 



w(S) 



tr c e 



-0(H(S)-nN) 



(5) 



Equation (0J is the starting point of Monte Carlo simu- 
lations of the Kondo models where the sum over the clas- 
sical spins is performed via importance sampling. The 
spin configurations S enter the Markov chain according 
to the weight factor w(S) that is computed via exact 
diagonalization of the corresponding one-particle Hamil- 
tonian in Eq. In the ID case we have performed MC 
simulations in which spins in domains of random lengths 
were rotated. We have performed MC runs with 2000 
measurements. The skip between subsequent measure- 
ments was chosen to be some hundreds of lattice sweeps 
reducing autocorrelations to a negligible level. 

As previously shownS^, the spin-integrated one-particle 
Green's function can be written as 



« o^jot, ».= / V[S] w(S)u}j(S) « c.;cj , 

a J 

(6) 

where <C c i ; ct is the Green's function in local spin 
quantization. It can be expressed in terms of the one- 
particle eigenvalues and the corresponding eigenvec- 
tors |-i/>( A ') of the Hamiltonian H(S): 



«c,;c] >£= 



It should be pointed out that the one-particle density of 
states (DOS) is independent of the choice of the spin- 
quantization. 



C. Uniform Hopping Approach 

The integral over the corespin states in the partition 
function |0J can be evaluated approximatively by resort- 
ing to a uniform hopping approach (UHA) 28 . The key 
idea is to replace the impact of the locally fluctuating 
corespins on the hopping amplitudes by some global av- 
erage quantity u. Then the Hamiltonian merely depends 




FIG. 1: Idealized spin and hole-density configuration in a ID 
Kondo chain at the critical chemical potential n* . A FM 
domain of Lf —4 lattice sites is embedded in an AFM back- 
ground. A single hole is localized in the FM domain giving 
rise to the depicted hole density (different from the schematic 
shape in Fig 4 of Ref.l30l). 



on one parameter, namely u, and the partition function 
can be written as the onc-dimensional integral 



Z 



duT(u) tv c e~^ A ^^ . 



The density of corespin states T(u) can be obtained nu- 
merically for 3D systems and analytically for ID systems. 

Obviously, this simplification assumes a uniform 
medium. In order to cope with magnetic polarons, the 
UHA has to be generalized and T(u) should be replaced 
by a two-parameter density T(uf,u a ), where Uf(u a ) de- 
notes the average hopping within the (anti)ferromagnetic 
domains. The details and results of such a generalization 
are the contents of Sec. lIVI of this paper. 



III. FERROMAGNETIC POLARONS 

Near half filling of a single e g band, a tendency towards 
phase separation has been observed in various studies. 
It has been claimed that the system separates into FM 
domains of high carrier concentration and AFM domains 
of low carrier concentration. In the following we show 
that a different picture rationalizes the ID Monte Carlo 
results in the range n w 0.7 — 1.0. 

We show that ferromagnetic polarons, i.e. single charge 
carriers surrounded by small ferromagnetic spin-clouds, 
are formed when holes are doped into the half filled 
e g band. In order to model such a polaron in a one- 
dimensional system, we take Lf adjacent lattices sites to 
be in ferromagnetic order and use Fl-l / (w) to account 
for the degeneracy of the remaining spins. 

First, we estimate the size Lf of the FM polaroid 
using a simple polaron picture. In this view the hole is 
confined in a perfectly FM domain consisting of Lf lattice 
sites and outside the domain the system is in perfect 
AFM order (see Fig.[T}. The tight-binding energies m a 
potential well (FM domain) with infinite barrier height 



are 



-2 cos i 



Lf + V 



(7) 
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The energy difference between a) a one-polaron state 
with perfect FM spins within the polaron and perfect 
AFM order outside and b) perfectly antiferromagneti- 
cally ordered ti g spins is given by 

Ae p = -2 cos (y~Tt) + 2 J <=ff ( L f - !) > 
\Lf + l/ 

where the first term accounts for the kinetic (derealiza- 
tion) energy of the hole in the potential well and the sec- 
ond term describes the energy deficiency due to (Lf — 1) 
ferromagnetic bonds. We have introduced the effective 
antiferromagnetic coupling J e g, given by 

Jeff = J' + 1/{2J H ) 

near n w 1 (see Ref. ITsI) . For typical values Jh = 6 
and J' = 0.02 we have J e ff ~ 0.1. Upon minimizing Ae p 
with respect to Lf, we obtain the optimal size of the 
polaron, which in the present case lies between Lt — 3 
and Lf — 4. If the FM domain contains N > 1 charge 
carriers, the energy difference is simply 

N 

Ae p (N) = -2 £ cos (/fy) + 2 Jeff (L/ - 1) . 
i/— l 

For N = 2, the optimum bi-polaron size is L/ ~ 7 and it 
increases to Lf ~ 10 for JV = 3 charge carriers. 

Next we estimate the chemical potential /j,*, at which 
holes start to populate the polaron states. Apart from 
the energy of the antiferromagnetic t2 g spins, the total 
energy (at T — 0) of the filled e g band is given in the 
grand canonical ensemble by —fi* L. By equating this 
energy to the total energy of the polaron we have 

-H*L = -fi* (L-1) + Ae p (l) 

which yields the desired chemical potential 

// = -Ae p (l) = 2cos (y^—) - 2J cff (L f - 1) . 
Lf + l 

The critical chemical potential is depicted in Fig. 
Similar considerations yield that /i* also approximately 
presents the limiting chemical potential between the filled 
antiferromagnetic band and a state with several single 
FM polarons, provided the polaron density is low, i.e. as 
long as we have an antiferromagnetic background. Con- 
sequently, at the chemical potential /i* the electron den- 
sity is not fixed. The energy gain Ae p exactly balances 
the loss of the chemical potenital — fj,*. This implies that 
the electron density has a discontinuity at fi* , i.e. the 
compressibility of the electrons diverges which has been 
previously interpreted as a consequence of PS tendencies 
(see Fig. 3a in Ref. |2(J and Ref. Il3|) . 

If we repeat the considerations for bi-polarons with the 
respective optimized size {Lf ~ 7 for the standard pa- 
rameter set), we find a critical chemical potential, also 
depicted in Fig. [21 which is very close to that of po- 
larons. If we proceed to tri-polarons we find again a 
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FIG. 2: Critical chemical potential \fj, c \ of the homogeneous 
model (solid line) from Ref. as compared to the chemical 
potential /i* at which polarons (dash dotted) and bipolarons 
(dashed) start to enter the Kondo chain. 

similar fx* . As can be seen in Fig. [21 the chemical po- 
tential fj,* virtually coincides with the chemical potential 
[i c of the "phase separation" obtained in Ref. 18. This 
potential was calculated as the point of coexistence of 
FM and AFM order in a homogeneous system. In view 
of the rather rough estimate of /x* we conclude that po- 
larons, bi-polarons, up to phase separated FM regions, 
are energetically comparable as long as the size of the FM 
domains is optimally adapted. For J e g > 0.12 which is 
comparable with the antiferromagnetic exchange in man- 
ganites, individual polarons are energetically favored. 

FM polarons, however, have a much higher entropy 
than the other objects and are therefore thermodynami- 
cally favored, even at low temperature. Moreover, at T ^ 
the domains are not completely (anti)ferromagnetically 
aligned which further reduces the energy differences 
between polaron and bi-polaron/phase-separated states 
considerably. Therefore, even for values of J c s < 0.12 we 
find individual polarons at very low temperatures. This 
conclusion is corroborated at (3 = 50 and J c s — 0.10 by 
the ensuing analysis of MC simulations. 

In order to scrutinize the polaron arguments, we com- 
pute the mean particle numbers for the corespin config- 
urations entering the Markov chain of unbiased MC 

tr (N e-PiHis)-^)) l 

(N) o •= tIc ^ v e L = V - 

V ' t r c (e-«^)-f*)) Aji + e/SMS)-/*) ' 

where e v (S) are the eigenvalues of H (S) for the con- 
figuration S. As a consequence of the above reasoning 
we expect a broad distribution of integer-valued par- 
ticle numbers if the chemical potential is close to /x*. 
The MC time series for (N)s f° r a L = 50 site chain 
(J H = 6, J' = 0.02, /3 = 50) is shown in Fig. EJ One time 
step corresponds to 1000 sweeps of the lattice. The in- 
verse temperature (3 = 50 corresponds to T ~ 50 — 100 K, 
i.e. a temperature relevant for experiments. 
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FIG. 3: Mean particle numbers in a grand canonical MC simu- 
lation (L = 50, J H = 6, J' = 0.02, P = 50) as a function of MC 
time. One time step corresponds to 1000 sweeps of the lattice, 
a) AFM case (fi = 1.22) b) polaron regime (/j, ~ /i* ~ 1.02) 
c) FM regime (/x = 0.80). 



The left-hand panel corresponds to a situation where 
the chemical potential is far above the critical chemical 
potential, which has the value p.* ~ 1.02 for the present 
parameter set. We see that the band is almost com- 
pletely filled, with isolated dots at N e = 49 correspond- 
ing to occasional FM polarons. At fi* (central panel), 
in agreement with the polaron picture we find a broad 
distribution of integer-valued mean particle numbers. If 
the chemical potential is reduced below fj,* , the system 
becomes ferromagnetic and we find the standard result 
of free electrons with a narrow and continuous spread in 
(N) not restricted to integer values. 

Hole-dressed spin-spin correlations provide another 
piece of evidence in favor of FM polarons. The bulk of 
Monte Carlo snapshots (not shown), as taken from sim- 
ulations for the FM Kondo model, contains isolated FM 
polarons of size Lf = 3 or Lf = 4. Once in a while two of 
them collide and form passing bi-polarons. The observed 
fraction of bi-polarons corresponds to a random distri- 
bution of polarons. In order to quantify the information 
revealed by the MC-snapshots, we introduce a modified 
corespin correlation function 



S n (l) 



1 



L-l 



L-l 



i+J 



(8) 



that measures the corespin correlations in the vicinity of 
a charge carrier (hole) . The density operator n\ for holes 
at site i is related to the density operator for electrons 
via = 1 — Tij. Figure 21 shows the results of an un- 
biased grand canonical MC simulation. The observables 
are evaluated at different subspaces with a fixed particle 
number. We observe ferromagnetic correlations that van- 
ish at / = 3 corresponding to a polaron that extends over 
Lf = 4 lattice sites. It should be pointed out that the 
MC result is almost independent of the number of holes 
in the system. In particular, the data do not indicate 
any enlargement of the FM domain for a larger number 
of holes. This result can only be explained by individual 
FM polarons because the size of the FM domain would 
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FIG. 4: Modified spin-spin correlation function from un- 
biased MC for an L = 50-site chain containing one (x), 
two (o), three (o), four (<), and five (>) holes. The inset 
shows the conventional spin-spin correlation function S(l) = 
j^j Y^,i=i Si ■ Si+i. The dashed line is calculated within the 
simple polaron picture, while the solid line represents the 
generalized UHA result for a single polaron. Parameters are 
P = 50, J' = 0.02, and J H = 6. 



strongly increase if there were two or more holes trapped 
in it. 

The inset of Fig. 0] shows the conventional corespin 
correlation function S(l). We observe the expected anti- 
ferromagnetic correlations, which decrease slightly with 
increasing number of holes. 

The result for the modified spin-spin correlation func- 
tion can again be explained qualitatively by the simple 
polaron picture. We consider a single polaron in which 
one charge-carrier is confined. Let the probability for 
the hole to be at site v in the FM region be p„, which 
is roughly given by the result for a particle in an infi- 
nite potential well (see Fig.^): p v cx sin 2 (^ ir/(Lf + 1)). 
The spin correlation is computed assuming perfect FM 
order inside the polaron and perfect AFM order outside. 
The result of this simple idea for Lf — 4 is shown as the 
dashed line in Fig. It agrees qualitatively with the MC 
data. 

For a quantitative, but still fairly simple description, 
we will now generalize the uniform hopping approxima- 
tion, to allow for FM polarons. 



IV. UHA FOR FM POLARONS 

In the previous section we have interpreted the MC 
data by using the simplest polaron ideas. In what follows 
we will refine our polaron picture by including thermo- 
dynamic fluctuations of the corespins. This is done by a 
generalization of the finite-temperature uniform hopping 
approach (UHA) introduced in Ref I281 

In the spirit of UHA, the impact of the corespins on 
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the motion of the e g electrons is now described by two 
UHA parameters, ut and u a , for the FM and AFM re- 
gion, respectively. These two parameters are averages of 
the hopping amplitudes in the FM and AFM domains, re- 
spectively. Their distribution is given by a two-parameter 
density of states denoted by TN f ,N a {uf, u a ). 

The size of individual polarons is fixed to Lf lattice 



sites. It is, however, possible that polarons overlap. The 
positions of the polarons are specified by the locations 
{ii, . . . ,i m } of their left ends, where m is the number 
of FM polaron wells. If polarons overlap, they may 
form a bi-polaron or even greater accumulations of holes. 
The grand canonical partition function in this generalized 
UHA reads 



Z = E E 



m /.. } 



du f du a T Nf , Na (uf, Ua ) tTce -0(Z(uf^>->^)-^) 



(9) 



r 



For ID chains, subject to open boundary condition, the 
joint density Tjf f ^ a (uf,u a ) depends merely upon the 
number Nt (N a ) of bonds in FM (AFM) regions. In 



higher dimensions, TN f ^N a (uf ,u a ) would actually depend 
upon the location of the individual polarons. The gener- 
alized UHA Hamiltonian reads 



H(uf,u a ; h, . . . ,i m ) = - ^ u ij 4 C 

<ij> 



Z 

1 ~2Jh 



r 



i+S ™i 



j' 4-i) 



(10) 



where 5 stands for n.n. vectors and z denotes the coor- 
dination number. The UHA parameter Uy are either Uf 
or u Q , depending upon the type of magnetic order at the 
adjacent sites i and j. The integrand of the partition 
function in Eq. © defines the joint thermal probabil- 
ity density p(uf, u a \ (3). From p(uf, u a \(3), we estimate 
mean values Uf,u a of Uf and u a of a L = 20-site chain 
reading 



0.937 , 



0.31 



(11) 



for the standard parameter set J' = 0.02 and Jh = 6 
and (3 = 50. These mean values are independent of the 
number of polarons and their positions as long as the 
total volume of the polarons is small compared to the 
system size. In order to simplify the following discussion, 
we will replace the thermal averaging of an observable by 
the value of that observable at v,f,u a - For the partition 
function this yields the simple form 



Z = 



E E 



tr,e 



- (3{H{uf ,u a ',ii , ■ • .,%)- /i(L— m)) 



(12) 

where the influence of the thermal fluctuations of the 
corespins is contained in the average hopping amplitudes 
Uf and u a - 

These average hopping amplitudes, however, are not 
sufficient for the determination of observables that do not 
directly derive from the partition function such as spin- 
spin correlations and one-particle spectral functions. In 



principle, these observables can be calculated in UHA by 
averaging over a set of typical thermal corespin configu- 
rations {S} obtained in UHA. In order to construct such 
a set, the azimuthal angles are also required, although 
they do not enter the energy and have a flat thermal 
probability density. The simplest way of constructing 
typical corespin configurations is to draw azimuthal an- 
gles at random. Starting from the reference spin S„ we 
proceed to the neighboring corespins by adding a random 
azimuthal angle x to the fixed relative polar angle. Thus 
we obtain a collection of typical corespin configurations 
{S}. 



A. Static Correlations 

We continue the discussion of the modified spin-spin 
correlation function S n (l). For a quantitative, but still 
fairly simple description, we take the deviations from per- 
fect FM and AFM order into account, while for the hole 
the approximate probabilities p v are retained. We em- 
ploy the mean UHA parameters Tif and u a to describe 
the relative angles of neighboring corespins and average 
over typical corespin configurations with random hole po- 
sitions. A comparison of the unbiased MC results with 
those of this approximate polaron approach, depicted in 
Fig. reveals an excellent agreement. 

Another observable to distinguish between polarons, 
bi-polarons, or even phase separated high density FM 
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clusters is the density-density correlation function 



L-l 

C V) : =T~iT, « n i - ( n i))( n li - • ( 13 ) 

i=l 



If holes form independent FM polarons, the correlation 
function should be structureless, while if holes gather in 
one FM regime, the correlation function will exhibit a 
positive peak at a typical inter-particle distance. In Fig.0 
C(l) is shown for the two-hole subspace, where only those 
spin configurations of the Markov chain are taken into 
account, for which Nh — 2. The UHA-polaron result is 
derived as follows. The positions of two FM-regions of 
size Lf = 4 are chosen at random, including overlap- 
ping ones. The hopping parameters are ut (u a ) for FM 
(AFM) bonds and the resulting tight-binding model is 
solved. The lowest two eigenstates will then be localized 
in the two FM potential wells. The resulting correlation 
functions are averaged over all possible positions of the 
FM potential wells. We observe a strikingly close agree- 
ment with the unbiased MC results. Similarly we proceed 
in the bi-polaron case, which is characterized by a single 
FM region of optimized size [L f = 7). Here the two holes 
occupy the ground state and first excited state of the FM 
potential well. The resulting correlation function differs 
drastically from the MC data. 
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FIG. 5: Density density correlation function for L = 50, /3 = 
50, J' = 0.01, Jh = 6, and Nh = 2. Error bars represent 
unbiased MC data, crosses stand for polaron-, and circles for 
bi-polaron results in UHA. 



This discrepancy increases with increasing hole num- 
ber, which shows clearly that the physics of the ID FM 
Kondo model is correctly described by single-hole po- 
larons and not by phase separation. 



B. Polaronic features in the spectral density 

In this subsection, we compute the spin-integrated 
spectral density 

Mu) = -- 5 « «L »- 

a 

for the original DE model by unbiased MC simulations 
based on the expression Eq. JJ]J. The results will again be 
analyzed by UHA in the framework of the polaron ansatz. 
In particular the pseudogap in the spectral density near 
n = 1, found in the FM Kondo modeli£*2ii, can readily 
be explained in the polaron picture. It is a consequence 
of the ferromagnetic box in which the hole moves. 

The computation of Ak(oj) in UHA is based on the 
reasoning that led to the partition function in Eq. I|12(l . 
I.e., for each polaron configuration and each value of u a , 
eigenvalues e' A ' and eigenvectors \ifj^) of the respective 
tight-binding Hamiltonian are determined, from which 
the Green's function 

is determined in local spin-quantization. The transforma- 
tion to the global spin-quantization is given by Eq. © 

er 

In the framework of UHA the relative angles of neighbor- 
ing spins are fixed by the parameters u a and Uf. For a 
unique description of the entire spin configuration, how- 
ever, azimuthal angles are again required. We proceed 
like in the discussion corresponding to Eq. (|13|l . i.e. for 
fixed parameters (uf,u a ), spin configurations are gener- 
ated with a flat sampling distribution in azimuthal an- 
gles. This is little numerical effort, as only the prefactor 
ujj (S) is affected. 

In the transition region (/x ~ fi*), the number of po- 
laron wells is not well defined, as pointed out in con- 
junction with Fig. In order to obtain a detailed un- 
derstanding we compare unbiased MC data and results 
of the polaron ansatz in the subspace of fixed particle 
(polaron) number. 

1. Antiferromagnetism at half filling 

We begin the discussion with the spectrum of the AFM 
state for the completely filled lower Kondo band (no po- 
larons). The unbiased MC data as shown in Fig. HO dis- 
play broad structures due to incoherent motion of charge 
carriers in a spin background that exhibits random de- 
viations from perfect AFM order. The result can be de- 
scribed by the UHA ansatz. For the mean value of the 
hopping amplitude we have u a — 0.31 which corresponds 
to (SiSi+i) = cos(#) ~ —0.7. Since there are no holes, 



FIG. 6: Spectral density of a half-filled Kondo chain with pa- 
rameters L = 20,(3 = 50, J' = 0.02, J H = 6, and u = 1.25. 
Dashed lines represent unbiased MC data, solid lines stand for 
UHA result (u a — 0.31). Inset: density of states. The verti- 
cal bar indicates the chemical potential. The left (right) -hand 
panel shows UHA-results in local (global) spin quantization. 
Error bars have been omitted for clarity. 



and consequently no polarons, the UHA parameters are 
the same for all bonds. The left-hand panel in Fig. 
shows the UHA result in local quantization, i.e. without 
the transformation given in Eq. (|14f) . In local quantiza- 
tion, the spectral density is simply given by 

Ak(w) = S(lu — 2u a cos(fc)) . 

The agreement with the unbiased MC result is rather 
poor at this stage although the band width is already well 
approximated. In the MC spectra, only very weak rem- 
nants of the tight-binding features are visible on top of 
the incoherent spectrum, which is almost k- independent. 

However, if we take into account the necessary trans- 
formation Eq. (|14|) to a global spin quantization, the 
agreement is strikingly close (see right-hand panel). The 
quasi-particle dispersion is strongly smeared out due to 
the random azimuthal angles of the corespins. There 
are two minor discrepancies between UHA and unbiased 
MC. The tight-binding remnants are more pronounced 
in UHA, while the MC results exhibit a weak structure 
near ui = 0, which is due to random fluctuations in the 
relative n.n. angles, resulting in locally trapped electrons. 
Nonetheless, at half-filling UHA and unbiased MC sim- 
ulations yield compatible results for the spectral density 
and the density of states (insets of Fig. |f)J| . 



2. Ferromagnetic Polarons 

Next we consider the case of one hole in the other- 
wise half-filled Kondo chain. To this end, we investigate 
the grand canonical MC data in the Nh — 1 subspace. 
The respective MC spectrum is shown in Figs. [7] and [S] 
The main feature of the spectrum is a broad incoherent 
background, similar to the one found in the AFM case. 
In addition, two dispersionless structures are visible at 
u> « ±1.5. As discussed earlie»i2i^S, a pseudogap shows 
up at the chemical potential. We find that an additional 
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FIG. 7: Spectral density for Nh = 1 holes (one polaron). 
Parameters like in Fig. [(J except /i = 1.1. Comparison of MC 
data (gray scale plot) with results of the simple polaron-well 
model. Circles (solid bright lines) stem from states localized 
in the FM (AFM)-well. Hopping parameters for the polaron- 
well model are Uf = 0.937, u a = 0.31, and Lf — 4. 



(mirror) gap appears at the opposite side of the spec- 
trum. 

Before discussing the UHA result, we want to provide 
a rough explanation of the MC data in terms of a simple 
polaron model. We assume an FM polaron-well of size 
Lf = 4, characterized by a tight-binding hopping param- 
eter Uf embedded in an antiferromagnetic background 
with hopping parameter u a . Since the two hopping pa- 
rameters are very different, we treat the various regions 
as separate chains and neglect their interaction. I.e. there 
is an isolated tight- binding chain of size L f corresponding 
to the FM region, and one or two chains corresponding 
to the AFM background. A carrier is localized either in 
the FM or the AFM domain. The eigenvalues in the FM 
region, e[ = — 2v,fCOs(kl), depend on the momentum 
kl = vixjiLf + 1) with v — 1, . . . ,Lf. The energies of 
states corresponding to different localizations of the FM 
well (io = 1, • • • ,L— Lf) are degenerate. The correspond- 
ing eigenstates in real space read 



tpl (X 



sin((i 
0, 



i + 



i < i < i + Lf 
otherwise. 



In this simple polaron-well model the expression for that 
part of the spectral density in local spin-quantization is 



L f L—L i 



y—l io — 1 



(15) 



£ Cv{k)8(u>-4) 



where \k) stands for the eigenvectors of the homogeneous 
tight-binding model with open boundary condition, i.e. 
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FIG. 8: Spectral density for Nh = 1 hole (one polaron). Sym- 
bols like in Fig. [(J and parameters like in Fig. |7| The broad 
incoherent central part of the spectrum (width ~ 2u a ) de- 
rives from the motion of the electrons in the fluctuating AFM 
background. Polaronic peaks show up at w ~ ±1.5. In be- 
tween these two structures a pseudogap opens at the chemical 
potential /i*. It is accompanied by a mirror pseudogap near 
-H*. 



(j\k) oc sin(jfc). The coefficient c u (k) as function of k 
shows a broad hump at k ~ fc„. Hence, the contri- 
butions of the FM regions to Ak(uj) are dispersionless 
structures at energies ui = el which are concentrated 
about k = k[. These structures are marked by open 
circles in the grayscale plot of Fig. [7\ They explain the 
additional features at the band edges, which are clearly 
visible in the spectral density in addition to the broad 
incoherent background. The latter is due to the mo- 
tion of the hole in the AFM regions. Since the AFM 
regions are much larger than the FM well, a continuous 
tight-binding band develops, characterized by the hop- 
ping parameter u a . The band is shown as a white line 
in Fig. [7| This part of the spectrum is similar to that at 
half filling (Fig. [fjj) . The remaining discrepancy as com- 
pared to the MC data is due to the fluctuations of the 
azimuthal angles of the corespins causing the white line 
to become more incoherent, as shown by the UHA cal- 
culations below. The transformation from local to global 
spin-quantization has, however, negligible impact on the 
polaron states in the spectrum, since they are due to the 
FM region, in which the fluctuations of the corespins are 
less relevant. As we can see in Fig.[7| our reasoning based 
on a single polaron well already describes the qualitative 
features correctly. 

The origin of the pseudogap and its 'mirror image' 
on the opposite side of the spectrum can now be sim- 
ply identified as the energy difference between the up- 



permost (lowest) state in the FM potential well [E = 
=F2u/cos(7r/(L/ + 1))) and the upper (lower) edge of the 
tight-binding band in the AFM region (E = ^2u a ) lead- 
ing to a width of the pseudogap AE = 2(v,f cos(7r/ (Lf + 
l))-u„), see Fig. □ 

A different picture would emerge if the FM domains 
were more extended, as for example in a PS scenario. 
They would then contain many energy levels (not only 
four as in Fig. CI) an d should thus not give rise to a pseu- 
dogap. 

For a more quantitative description we invoke the two- 
parameter UHA as described before. In Fig. [S] the UHA 
results, already in global quantization, are compared with 
those of unbiased MC simulations. The features of the 
spectral density are well reproduced. UHA even yields 
quantitative agreement as far as the pseudogap in the 
density of state is concerned. 




FIG. 9: Spectral density for Nh = 2 holes (two polarons) in 
the left-hand panel and Nh ~ 3 holes (three polarons) in the 
right-hand panel. Parameters and meaning of symbols like in 
Fig. H 

For the same value of the chemical potential (i = 1.1, 
the two- and three- hole subspaces have been studied by 
projecting the grand canonical MC data on the Nh — 
2 and Nh = 3 subspace, respectively. The results are 
depicted in Fig. |5J The spectra are qualitatively similar 
to those for the one-hole subspace. Only the spectral 
weight of the polaron peak increases and shoulders show 
up in the AFM part of the spectrum. They are due to 
the interstitial AFM regions enclosed by the FM-polaron 
wells. The argument is the same as before. The allowed 
energies are e(k) = — 2u a cos(fc), but now the possible 
fc-values depend on the size of the interstitial regions. 

C. Discussion 

The emerging global picture is as follows. There ex- 
ists a critical chemical potential ji* . The value of fj,* can 
be obtained from simplified energy considerations. For /i 
significantly above fi*, the band is completely filled and 
antifcrromagnetic. The spectral density in this case is 
shown in Fig. |SJ At fi*, according to Fig. El holes enter 
the e g -band forming isolated FM domains each contain- 
ing a single hole (see Figs. 0] and [SJ|. In the grand canoni- 
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cal ensemble, the number of polarons strongly fluctuates 
and the height of the polaron peak in the spectrum is 
directly linked to the number of holes. The pseudogap 
appears around the critical chemical potential. For val- 
ues below ji* , the system switches from predominantly 
AFM order to FM behavior and the pseudogap gradu- 
ally disappears in favor of a single quasi particle band 
of tight-binding type. Our analysis yields compelling ev- 
idence against the PS scenario and in favor of FM po- 
larons in ID. Furthermore, it appears plausible that the 
formation of FM polarons will exist in symbiosis with 
lattice deformations (Jahn- Teller polarons )£i2£. 

For CMR oxides, Saitoh et al— have investigated the 
temperature dependence of angle resolved photoemission 
spectra (ARPES) for the phase transition from FM to 
PM order. These studies show that a pseudogap also 
develops above Tc, which can be rationalized in the po- 
laron picture. In the PM phase we have the competi- 
tion of ferromagnetism, driven by the DE mechanism, 
and spin disorder due to thermal fluctuations. There- 
fore, FM polarons will form in the paramagnetic back- 
ground. They will, however, be more extended because 
the PM force is less pronounced than the AFM force 
at low temperatures^. With increasing temperature, the 
corespin fluctuations become stronger and the competi- 
tion of the FM polarons with the PM background gets 
tougher. The existence of FM domains above Tc has 
been corroborated by neutron scattering experiments^. 

Furthermore, the ARPES experiments revealed that 
the bandwidth changes merely by about 4% across the 
FM to PM phase transition. On the other hand, it has 
been argued^ that the DE model predicts a reduction 
of about 30% in the PM phase if there the mean angle 
between neighboring spins is taken to be 7r/4 and the 
mean hopping parameter is therefore reduced to 1/V2- 
The authors in Ref. |53 therefore conclude that 'DE is 
probably not even the dominant mechanism . . . '. 

At first glance, the argument seems convincing. How- 
ever, in the polaron picture we do not really expect such a 
dramatic change of the band width since it is determined 
by the polaronic peaks at E = ±2u f cos(ir / (L f + 1)). 
Hence, the band edges depend on the hopping parameter 
of the FM region and not on that of the PM region. More- 
over, by the same reasoning that leads to a bandwidth 
reduction of 30% in the PM state, one would expect that 
the bandwidth vanishes in the AFM phase at low temper- 
atures, since here the neighboring corespins are mostly 
antiparallel. That conclusion is in strong contrast to the 
unbiased MC results depicted in Figs. |S1 and [SJ Even for 
the incoherent inner part of the spectrum, which is due 
to electronic motion in the AFM region, a considerable 
band width exists, due to the spin fluctuations which 
are present even at very low temperatures. On top of 
that, the band-edges at finite hole filling are not really 



determined by the AFM regime, but rather by the FM 
polarons. 

V. CONCLUSIONS 

In this paper, polaronic aspects of the ferromagnetic 
Kondo (double-exchange) model have been analyzed by 
unbiased finite temperature Monte-Carlo simulations and 
they have been explained by simple physical pictures. It 
has been found that in ID, the physical effects of the 
FM Kondo model close to half filling are not governed 
by phase separation, as previously reported, but rather 
by single-hole ferromagnetic polarons. They can be ex- 
plained qualitatively on the back of an envelope by ide- 
alized polaron pictures. 

It seems sensible to reassess the explanations of CMR 
based on PS. These explanations are primarily based on 
percolation ideas, which can equally well be applied to 
FM polarons as percolating units. It appears plausible 
that the formation of FM polarons will exist in symbiosis 
with lattice deformations (Jahn- Teller polarons). Single- 
hole FM polarons allow a direct explanation of the pseu- 
dogap, observed in the manganites, whereas for larger 
FM clusters the pseudogap would be filled up by addi- 
tional states. The striking similarity of the bandwidth 
of the FM and the PM phase, observed in ARPES ex- 
periments, can also be explained by FM polarons in the 
frame of the DE model. Moreover, the infinite compress- 
ibility near the half filled band, which has previously been 
attributed to PS, is a consequence of the fluctuating num- 
ber of polarons in the grand canonical ensemble. Work 
is in progress for higher-dimensional systems where the 
entropy is expected to have less influence on the thermo- 
dynamic behavior. 

For the analysis of the Monte Carlo results, we have 
extended the uniform hopping approach (UHA) at finite 
temperatures to include polaronic effects. This ansatz 
reduces the numerical effort by several orders of magni- 
tude, while retaining all crucial physical features. The 
key idea is to map the physics of the high-dimensional 
configuration space of the ti g corespins onto an effective 
two-parametric model. A full thermodynamic evaluation 
of the UHA model takes into account entropy and fluctu- 
ations of the corespins. The results are in close agreement 
with the unbiased MC data and allow a realistic descrip- 
tion of all FM polaron effects found in various physical 
quantities. 
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